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Abstract. We present optical to very-high energy 
(VHE) gamma-ray observations of Mrk 421 between 
2008 May 24 and June 23. A high-energy (HE) 
gamma-ray signal was detected by AGILE-GRID 
during June 9-15, brighter than the average flux 
observed by EGRET in Mrk 421 by a factor of 
« 1.5. In 20-60 keV X-rays, a large-ampUtude 5- 
day flare (June 9-15) was resolved with a maximum 
flux of approx. 55 mCrab. Super AGILE, RXTE/ASM 
and Swift/BAT data show a clearly correlated flaring 
structure between soft and hard X-rays, with a high 
flux/amplitude variabUity in hard X-rays. Hints of 
the same flaring behavior is also detected in the 
simultaneously recorded GASP-WEBT optical data. 
A target of opportunity observation by Swift near the 
flare maximum on June 12-13 revealed the highest 
2-10 keV flux ever observed (>100 mCrab) and a 
peak synchrotron energy of approx. 3 keV, a large 
shift from typical values of 0.5-1 keV. Observations 
at VHE (E >200 GeV) gamma-rays during June 6-8 
show the source flux peaking in a bright state, weU 
correlated with the simultaneous peak in the X-rays. 
The gamma-ray flare can be interpreted within the 
framework of the Synchrotron Self Compton model 
in terms of a rapid acceleration of leptons in the jet. 
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I. Introduction 

Mrk 421 is a nearby blazar {z = 0.031), detected in 
7-rays by EGRET |1| and it was the first extragalactic 
object detected at E > 500 GeV |2|. It belongs to the 
class of high-energy peaked BL Lac objects (HBLs) |3|, 
i. e. radio-loud active galactic nuclei with high radio 
and optical polarization. Its spectral energy distribution 
(SED) is double-humped with a first peak usually in 
the soft to medium X-ray range, and a second one at 
GeV- TeV energies |i4J, [SJ. The first hump is commonly 
interpreted as due to synchrotron radiation from high- 
energy electrons in a relativistic jet, while the origin of 
the second peak is still uncertain. In leptonic scenarios 
it is interpreted as inverse Compton (IC) scattering 
of the synchrotron (Synchrotron self-Compton, SSC) 



or external photons (External Compton, EC) by the 
same population of relativistic electrons. The observed 
correlated variability between X-rays and TeV 7-rays 
B, Q, 10 is well explained in the SSC framework 
whereas the EC scenario is unlikely to apply in HBLs, 
due to the low density of ambient photons. Alternatively, 
hadronic models invoke proton-initiated cascades and/or 
proton-synchrotron emission IfTOl . ifTTI . Leptonic and 
hadronic scenarios for HBLs predict different properties 
of the 7-ray emission in relation to emissions in other 
energy bands. 7-ray observations of flaring BL Lacs and 
simultaneous multiwavelength data are thus the keys to 
investigating these two scenarios. 

A hard X-ray flare of Mrk 421 was detected by 
SuperAGILE on 2008 June 10 fill. This detection 
was later followed by a detection in 7-rays ifTsll by 
the AGILE/GRID (Gamma Ray Imaging Detector) and 
prompted a ToO observation by Swift/XRT, comple- 
menting the ongoing multifrequency observing cam- 
paign of Mrk 421 with GASP-WEBT (optical), MAGIC 
and VERITAS (TeV). A full account of this work has 
recently been published [14|. 

II. Observations and Results 

A. AGILE observations 

The AGILE 1151 composite payload allows for si- 
multaneous observations in the energy ranges 30 MeV- 
30 GeV and 20-60 keV over a very wide field of 
view by means of GRID and the hard X-ray imager 
SuperAGILE, respectively. Mrk 421 was observed for 
five days, between 2008-06-09 17:02 UT and 2008-06- 
15 02:17 UT 

1) Hard X-ray observations: On 2008 June 10, Su- 
perAGILE detected enhanced hard X-ray emission from 
Mrk 42 1 . The measured flux in 20-60 keV was found to 
be above 30 mCrab (24-hour average), almost an order 
of magnitude larger than its typical flux in quiescence. 
In the following days, the flux increased up to about 
55 mCrab. The 5-day 20-60 keV SuperAGILE light 
curve is shown in Fig. Ic. Using the publicly available 
light curves for this source from the BAT instrument 
in the 15-50 keV energy range, we calculated daily 
averages, obtaining a good coverage also before and 
after the AGILE observation (black squares in Fig. 
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Ic), revealing that SuperAGILE indeed observed the 
maximum brightness of this hard X-ray flare. 

SuperAGILE photon-by-photon data allows extraction 
of a time-averaged energy spectrum from the mask- 
convolved data. Given the lack of substantial spectral 
variability in the SuperAGILE/ASM hardness ratio (Fig. 
le) we accumulated the average energy spectrum from 
the data of the last 4 days of the observations, when the 
source flux varied between 35 and 55 mCrab, for a total 
exposure of 140 ks, constraining the photon index of 
a simple power law, T = 2A3tosl (xLf = 0-8. 2 dof). 
The average flux is F(20 - 60keV) = (4.90 ± 0.54) x 
10-1° erg cm-2s-i ((9.8 ± 1.1) x lO"^ photons cm^^ 
s-i). 

2) Gamma-ray Observations: Mrk 421 was not de- 
tected on daily time scales, implying a daily-averaged 
flux < 100 X 10^^ photons cm^^ s^^, similar to what 
observed by EGRET |16|. A 4.5-(t significance in the 
range 100 MeV-10 GeV resulted from an integration 
over the whole 5-day period (^ 260 ks). The measured 
flux is (421:^2) X 10"* photons cm^^ g-i^ ^^0^ ^3 
times higher than the average flux detected by EGRET 
(r^ 13 X 10"®, fT6l) and --1.5 times higher than, but stifl 
consistent, with the highest flux (27±7) x 10^^ photons 
cm^-^ s^i observed during the Viewing Period 326. 

B. The Soft X-ray band 

1 ) Swift X-ray telescope: Following the SuperAGILE 
detection, on 2008 June 12 Swift ToO observations were 
triggered, between 2008 June 12 19:33:20 UT and June 
13 at 01:57:37 UT for 5 ks. The XRT spectral data are 
well described by an absorbed log-parabolic model. A 
joint fit of the XRT and SuperAGILE spectral data using 
the 4-day average spectrum used the log-parabolic model 
with Galactic absorption (N^'^^ = 1.61 x 10^° cm^^; 
iflTll ). defined as: 

F{E) = /s:£;-«-''iog(-E) photons cm^^ g-i 

where a is the photon index at 1 keV and b is the 
curvature parameter ifTSl . lfT9l . This model usually 
describes adequately the featureless and curved spec- 
trum in HBLs. The joint fit provides a = 1.65^0 q2' 
b = 0.37tom5 (xIof=l-4, 763 dof), which impUes a 
peak energy 2.97iQ'^g keV, and predicts F2_iokcV = 
2.56 X 10"^ ergcm"^ s-i (or 0.4 photons cm^^^ s^^) 
and F2o-60kov ^ (5.7 ± 0.6) x 10^1° erg cm^^ g-i 
((1.1 ± 0.1) X 10^^ photons cm^^ s^i), comparable to 
the stand-alone SuperAGILE best fit. 

2) RossiXTE All Sky Monitor (ASM): We retrieved 
the public light curves provided by the ASM to trace 
the evolution of the soft X-rays during the AGILE 
observation. Fig. lb shows the daily light curve of Mrk 
421 in the energy range 2-12 keV, obtained by properly 
averaging the dwell-by-dwell data. 

The emission at soft X-rays is well correlated with 
hard X-ray emission. The ASM data show that the 
XRT observation took place when the source was at 
its maximum emission at soft X-rays (MJD ~ 54630). 



Comparing the relative intensity of the two flares in Fig. 
1, the second peak appears to be significantly harder 
than the first one. This is also seen in Fig. le, where 
we computed the daily-averaged hardness ratio between 
hard (15-60 keV) and soft (2-12 keV) X-rays. The 
source appears to have undergone the hardest part of 
this double-humped flare just during the AGILE 7-ray 
detection. 

C. Observations in the VHE ^-rays 

The 4-telescope array VERITAS [20] in Arizona, 
USA, and the single-dish instrument MAGIC EH, ED 
at La Palma are imaging air Cherenkov telescopes cov- 
ering an energy range from ^100 GeV to some tens of 
TeV The instruments have a typical energy resolution of 
<20% (VERITAS) and 20-30% (MAGIC), and event- 
by-event angular resolution of < 0.14°. Wobble-mode 
observations ll23ll at an 0.4° offset from the camera 
center were taken on five nights from 2008 June 3- 
8 with MAGIC at zenith angles (ZA) between 28° 
and 48° and on 4 nights (May 27, June 1, 5, 6) wifli 
VERITAS (wobble offset: 0.5°) at ZA between 32° 
and 40° during partial moon light conditions. The total 
live-time after applying quality selection is 2.95 and 
1.17 hours with MAGIC and VERITAS, respectively. 
The data were analyzed using the MAGIC 1241 and 
VERITAS l25l, EOl standard calibration and analyses 
and image parameters ll26l . Two independent analyses 
of both the MAGIC and VERITAS data sets yielded 
consistent results. In total, a signal corresponding to a 
significance level of 44cr (VERITAS) and 66cr (MAGIC) 
is obtained. The combined MAGIC- VERITAS data (Fig. 
Id) show a transient peaking near MJD 54622. The 
VERITAS energy spectrum for June 6 is provided. A 
power law fit over the energy range 0.3-5 TeV resulted 
in a Xdof=0-7 with a photon index T = 2.78 ± 0.09. 
In Fig. 2 we show the intrinsic 7-ray spectrum at the 
source reconstructed by removing attenuation effects by 
the extragalactic background light II27I following the 
procedure of |28|. Fitting a power law to the intrinsic 
spectrum yields a photon index T = 2.59 ± 0.08, which 
is not significantly harder than the measured spectrum 
due to the relatively low redshift z — 0.031. 

D. Optical and UV observations 

Mrk 421 is one of the 28 7-ray-loud blazars regu- 
larly monitored by the GLAST- AGILE Support Program 
(GASP; L29J) of the Whole Earth Blazar Telescope 
(WEBT). The GASP observations started in early 2007 
November. 

During the Swift pointing on 2008 June 12-13, the 
UVOT Oni insti-ument observed Mrk 421 in the UVWl 
and UVW2 photometi'ic bands. The UVOTSOURCE 
tool is used to extract counts, correct for coincidence 
losses, apply background subtraction, and calculate the 
source flux. 
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III. Discussion 
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Fig. 1. a) i?,-band optical liglit curve from GASP-WEBT (May 24- 
lune 23); b) ASM (2-12 keV) light curve (bin size is 1 day) and XRT 
(2-10 keV) flux (triangle); c) SuperAGILE (20-60 keV, triangles; 1 
Crab = 0.2 ph cm-^ s-^) and BAT (15-50 keV, empty squares; 1 
Crab = 0.29 ph cm-2 s"!); d) MAGIC and VERITAS (>400 GeV, 
empty squares and filled circles, respectively), the Crab flux > 400 
GeV (horizontal dashed line), AGILE (> 100 MeV, triangle and right 
vertical axis); e) the hardness ratio computed by using the SuperAGILE 
and ASM data for each day. The shaded areas mark period 1 and period 
2. 



Mrk 421 showed an interesting broad-band activity 
during the first half of 2008 June as derived from AG- 
ILE data combined with those of GASP-WEBT, Swift, 
MAGIC and VERITAS. Using our multi-frequency data 
we were able to derive time-resolved SEDs (Fig. 2). We 
distinguish two time periods: period 1, 2008 June 6, 
and period 2, 2008 June 9-15. The optical, soft and 
hard X-ray bands strongly constrain the SED around 
the synchrotron peak, and its daily variability reveals 
the physical processes of Mrk 421. Possible correlated 
variability is shown in Fig. 1 between the optical, the 
X-ray, and the high-energy parts of the spectrum. Based 
on the physical constraints obtained for the synchrotron 
peak, we can model both the HE and VHE 7-ray emis- 
sion. We first model the synchrotron peak of emission 
using the period 1 optical, soft and hard X-ray data. 
The short time- variability (Fig. 1) constrains the size of 
the emitting region to i? < cT5 5 x IQ^'^{5/2Q) cm. 
Hence, we consider a one-zone SSC model |31| based 
on a blob of comoving size i? — 4 x 10^^ cm, with 
a relativistic Doppler factor 5 = 20 and characterized 
by non-thermal comoving electron energy distribution 
function described by a double power law: 

Ml) = (1) 



(7/76 



(7/7b)'' 



where the comoving Lorentz factor (7) varies in the 
range 7„.„i = 4 x lO^ < 7 < ^^.^x = 1.3 x 10^, the 
normalization (density) constant K — 4 x 10^^ cm^'^, 
and the break energy 7;, — 3.6 x 10^ and with pi — 2.22, 
P2 = 4.5, the low-energy and high-energy power-law 
indices, respectively (see Table 1). With these parameters 
we found that the data for period 1 are best fitted with 
a comoving magnetic field B — 0.1 G. 

Variability may be caused by several factors; we con- 
sider two cases: (A) hardening/softening of the electron 
energy distribution function caused by particle acceler- 
ation processes; (B) increase/decrease of the comoving 
particle density, as a consequence of additional particle 
injection/loss by shock processes. 

We expect TeV variability to be comparable with 
the X-ray one if case (A) applies: this is because the 
emission is in the Klein-Nishina regime. Alternately, for 
case (B) we expect the TeV relative variability (AF/F) 
to be a factor of 2 greater than that of the X-ray flux 
variability. 

Our AGILE, MAGIC and VERITAS data appear to 
support case (A). We compare the SEDs for period 1 
and period 2, to better assess the spectral evolution. In 
Fig. 2 we show our optimized modelling of the time- 
resolved synchrotron peak and consequent SSC high- 
energy emission for the period 1 as well as for 2008 June 
12-13 of the period 2. In the last case, the adopted model 
parameters are pi = 2.1, p2 = 5, 7b = 4.2 x 10*^, K — 
6 x 10^^ cm~'^. Our model predicts an even larger TeV 
flux for period 2 (no TeV observations exist, however) 
than detected in period 1. 
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Fig. 2. SEDs of Mrk 421 obtained by combining the GASP-WEBT, 
SWIFT/UVOT, Ros.siXTE/ASM, XRT, SuperAGILE, BAT, GRID and 
VERITAS data in period I and period 2 (empty and filled circles, 
respectively). Both are one-zone SSC models (dashed line for period 
I and solid line for period 2). 



The optical light curve shows variations of the order 
of 10% on a time scale '^few days, superimposed on 
a long decay during the entire period. Individual soft 
and hard X-ray peaks result in increased fluxes by a 
factor of ^^2.5 and ^^5, respectively: no long term decay 
appears. This different behavior of the X-ray radiation 
and the bulk of the optical emission may suggest a 
more complex scenarios than A) and B): optical and X 
radiation comes from two different jet regions, each one 
characterized by its own variability. A possible scenario 
is one in which the inner jet region would produce the 
X-rays and it would be at least partially transparent to 
the optical radiation. In contrast, the outer region can 
only produce lower-frequency emission. The signature 
of the X-ray events visible in the optical light curve 
would come from the inner region and would be diluted 
by the optical radiation emitted from the outer region, 
see ||32l for the case of Mrk 501; |[33l. 

Interestingly, the 2-10 keV flux measured by XRT on 
June 12 - 13, ^ 2.6 x 10"^ erg cm"^ s"^, is higher than 
all previous observations (< 2 x 10"^ erg cm"^ s"^; 
Q, in ). A joint analysis of the XRT and SuperAGILE 
data, covering the range from 0.7 to 60 keV, provides 
a best-fit spectral model consistent with a log-parabolic 
shape, with parameters implying a peak energy ^ 3 keV, 
in good agreement with the steeper positive correlation 
between the peak energy and the maximum of the SED 
found by f35|. although our value of the peak energy 
shows a significant shift with respect to typical values 
of 0.5-1 keV for this source. 
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